Can Binary Bound States in a Strongly Coupled Quark-Gluon Plasma 
be observed via dileptons and photons? 
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Recently there was a significant change of views on physical properties and underlying dynamics of 
Quark-Gluon Plasma at T = 170 — 350 Me V^, produced in heavy ion collisions at RHIC. Instead of 
weakly coupled gas of quasiparticles, it is rather a liquid-hke matter with multiple bound states. In 
this paper we discuss how one can test these ideas experimentally, using the "penetrating probes" 
and looking for certain peaks at some invariant masses. In dileptons the most promising are modified 
p,uj, with M{T « Tc) ~ .5 Gey and also near zero binding at M{T « (1.5 - 2)Tc) = 1.5 - 2Gey. 
We also discuss the observability of peaks corresponding to scalar/pseudoscalar mesons in the two- 
photon channel. 



I. INTRODUCTION 
A. The penetrating probes 

To use photons and dileptons as "penetrating probes" 
to study the earlier stages of heavy ion collisions was 
one of the first suggestions to study quark-gluon plasma 
(QGP) [1]. Although these particles are produced much 
less copiously compared to hadrons, and thus they are 
more difficult to study, photons and leptons are of special 
importance because they do not suffer final state interac- 
tion and therefore can bring us direct information about 
the earliest (the hottest) stages of the collision. In this 
respect one may compare them to e.g. solar neutrinos, 
which bring direct information about the conditions at 
Sun's interior. 

It was argued in [1] that because of the QCD phase 
transition, one should expect the invariant mass (M^ = 
(p;+ + pi-)'^) spectra of penetrating probes produced in 
QGP to be qualitatively different from those originating 
in hadronic matter. However only monotonously decreas- 
ing mass spectra were expected from QGP, in contrast 
to familiar vector meson p, w, 0, J/i/'... peaks from the 
hadronic phase. The rates were calculated in a picture of 
weakly interacting and freely moving near-massless quark 
quasiparticles, below wQGP for brevity. 

Recently a radically new picture of QGP has been de- 
veloped, to be described in next subsection, based on 
much stronger coupling between quasiparticles, to be re- 
ferred to as sQGP. One aspect of it is that the meson- 
like bound states continue to exist at T = (1 — 2)Tc = 
170 — 350 MeV (the temperature range corresponding to 
QGP at RHIC), although in a strongly modified form. 
Furthermore, even at higher T when there are no bound 
states, there exist near-threshold enhancement which 
may still be used to infer the value of the quasiparticle 
masses and the strength of their interaction. 

As penetrating probes are emitted during the whole ex- 
pansion of the fireball, the evaluation of observable spec- 
tra includes two very different steps, (i) First one eval- 
uates a production rate in equilibrium matter at a given 



T, dR{T) / d^xd^q per time per volume per 4-momentum. 
The second step (ii) is the space-time integration over the 
4-volume of the expanding fireball, from its creation to 
final freezeout. At the first step (i) there will be certain 
peaks at fixed values of mesons masses, but since they 
are in general T-dependent, my(r), the time integral at 
the second step will in general smoothen them out. Only 
peaks corresponding to special points can survive and be 
observable, as we detail below. 

Let us start with the basic expression for the dilepton 
production rate 
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where a is the electromagnetic coupling, q^ is the dilep- 
ton energy and the "formfactor" F for two well-tested 
processes are. 



dcf 
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F = 



where the former one corresponds to tttt — > p l^l^ 
annihilation in a vacuum or hadronic phase at small 
T, written in standard vector-dominance form. The 
latter expression corresponds to qq annihilation simi- 
lar to partonic Drell-Yan process. If one ignores quark 
masses in Fq, it is just a sum of squared charges for 
all relevant quarks, u,d,s of all colors. This basic case 
we will use as our "standard candle" below, normaliz- 
ing all predictions to the "standard wQGP rate" with 
PwQGP — 24. Detailed calculation of the final dilepton 
spectra for such "standard candle" model were done by 
Rapp and Shuryak [2] (including specific acceptance of 
the SPS NA50 experiment) and for RHIC by Rapp [4]. 
For a review see [5] . 

In terms of the imaginary part of the photon self energy 
wQGP limit means 
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If quark quasiparticles in the QGP phase arc not mass- 
less (and current lattice calculations indicate they are as 
heavy as IGeV in the region of interest Tc — 3Tc) one 
would like to correct for that using the formulae from 
above. This fact leads to a trivial modification to the 
previous expression (that is relevant in the IMR) 



9n= i 
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[ (M^ > 4m^) 



Furthermore, we expect nontrivial modification of the 
annihilation cross section of quarks into leptons due to 
the attractive interaction between them. This modifica- 
tion will be specially important for near threshold pro- 
duction, where we will use non relativistic Green's func- 
tions. Modifications of the rates are also expected due 
to the presence of bound states, where, as we sill see, the 
problem is intrinsically relativistic. 

Dileptons produced by vector resonances (p, iv. 0, J/ip) 
are described by the first expression in (2), which can 
be used both in a hadronic phase and in sQGP. There 
exist large literature on matter-induced modification of 
mesons in hadronic phase. For sufficiently dilute mat- 
ter one can use linear density approximation, in which 
the modification of real and imaginary part (the mass 
and the width) of a state is simply determined by for- 
ward scattering amplitudes in vacuum which are often 
experimentally known. For example, irp and Np scat- 
tering used in [7] are mostly responsible for modification 
of p resonance. Quite a number of works following such 
ideas found relatively modest shifts of resonance masses 
(e.g.TOp) downward and some broadening*. Further as- 
sumptions arc needed when the matter is not dilute: e.g. 
the Li-Ko-Brown model [6] assumed that all hadronic 
masses are shifted proportional to Walecka-type scalar 
mean field, related to local matter density. As exam- 
ple of self-consistent approach see e.g. a set of coupled 
gap equations discussed by Rapp and Wambach [5]: it 
predicted very strong broadening of the p meson due to 
mixing between p with excitations of the lowest baryon 
resonances, such as^ A^*(1520)iV~^. 

Let us now proceed to a brief review of dilepton exper- 
iments. Since these kind of measurements are generally 
much more difficult to do (there are large background to 
be rejected, while the process itself has small cross sec- 
tion), they have a sad tendency to come "too late", right 
before the program gets closed. 



At Berkeley BEVALAC the DLS dilepton spectrometer 
has found a spectacular effect, but the detector efficien- 
cies were debated years after its last run-t. Brookhaven 
AGS program had no dilepton experiments, while CERN 
SPS program included: (i) CERES (NA45), which stud- 
ied the low mass (M=0-1. GeV) e+e^ pairs, (ii) 
HELIOS-3, which studied medium mass /Lt+/x^" M=l-2 
GeV, and (iii) NA38/50/60 concentrated on p^p~ with 
M=2-4 GeV. All of them observed quite significant en- 
hancement over "standard hadronic sources", ranging 
from factor 5 at CERES (in some kinematic region) to 
about 3 at NA38/50 for intermediate masses § M=2-3 
GeV. The observed effective "p melting" at CERES and 
thermal continuum at NA50 indicates an approach to- 
ward chiral symmetry restoration and QGP [2,3]. This 
is quite puzzling, since at the SPS only a small fraction 
of the space-time volume contributing to dilepton pro- 
duction is expected to originate from QGP, while most 
of it should still be from the mixed or hadronic phase. 
The NA50 experiment now evolved into NA60 which is 
supposed to resolve the remaining uncertainty between 
thermal dileptons and charm contribution, as well as sig- 
nificantly increase mass resolution and tell us about pos- 
sible modifications of narrow oj, (j) resonances. 

The dileptons are supposed to be the main part of the 
RHIC program: PHENIX was aimed at that from the 
beginning and now also STAR detector has electromag- 
netic calorimeter. Large-statistics Run-4 data are now 
being analyzed, and the only published data so far are 
PHENIX single electron production [10], which is pre- 
sumably dominated by charm decays. 

Matter modification of more narrow ui, </> resonances 
can obviously be observed only for a fraction of these 
particles decaying inside the fireball. Resonances de- 
caying after freezeout have the usual vacuum mass, and 
these decays produce large "hadronic background" for 
such measurements: in fact so far none of the existing 
dilepton experiments had sufficient resolution to observe 
these modifications. However for sQGP vector mesons 
the mass modification is very large, e.g. we will discuss 
below peaks with masses up to M ~ 2 GeV, and so in 
this case the resolution is not so much an issue. 

In principle, the 77 channel of the decays provides an- 
other window into early-time spectroscopy, as photons 
are penetrating probes as well. Still, not much discussion 
of those can be found in literature. In-matter modifica- 
tion of the pseudoscalars tt, are expected to be quite 
small, as these masses are protected by chiral symmetry 
and are not expected to be modified much in all hadronic 
phase T < Tc- The experimental observability of these 



*As shown in [8], mass reduction leads to some suppression 
of the decay, which often tends to cancel widening from a 
rescattering. This was experimentally confirmed by STAR 
experiment at RHIC [9] , which found no width change of the 
p observed via nn at late stages. 

^This notation means a nucleon excited into a baryonic res- 
onance by vector meson absorption. 



* There are no ofBcial pulication yet, but preliminary results 
from HADES experiment at GSI recently reported have not 
confirmed this effect in its first run. 

^ At high masses region, M > 3 GeV, the dilepton production 
is well described by simple partonic Drell-Yan process. 
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shifts in 77 channel is next to impossible, dne to huge 
background from similar decays after freezeout. Scalars 
and tensors are considered to be too wide even in mat- 
ter, to get any meaningful signal even before in-matter 
modification. 

The situation again should bo quite different for 
7r°, T], rj' and scalars originating from sQGP**, since very 
large mass shifts are expected. We return to their dis- 
cussion in section IV below. 



B. Bound states in QGP and the penetrating probes 

The quark-gluon plasma (QGP) was defined as a state 
of matter in which a color charge is screened [1] rather 
than confined. According to lattice and experimental 
results it takes place above the critical temperature Tc w 
170 MeV. 

At high temperatures T 3> T), it should be a gas 
of weakly interacting quasiparticles (modulo long-range 
magnetism), wQGP. However a physical picture of QGP 
structure at not- too-high T, in the temperature range 
(1 — 3) Tc, recently underwent a radical change [11]. The 
interaction seems to be sufficiently strong to produce 
multiple bound states of quark and gluon quasiparticles, 
and wc will refer to such matter as strongly coupled QGP 
(sQGP for brevity). 

Using lattice-based potentials, it was shown [12,13] 
that a picture of binary bound states may provides a 
consistent description of several previously disconnected 
lattice observations, such as (i) bound states for char- 
monium and some light qq states [14]; (ii) static poten- 
tials [15]; (iii) quasiparticle masses [16]; and (iv) bulk 
thermodynamics [17]. However more complicated struc- 
ture can also be present, perhaps in line with a liquid-like 
behavior needed to understand robust collective flow phe- 
nomena at RHIC, well described by ideal hydrodynamics. 

Studies of strong coupling and related binary states are 
also under way in other fields of physics. In particular, 
the 7V=4 supersymmetric Yang-Mills theory at finite T 
is a perfect toy model, in which a parametric transition 
from weak to strong coupling can be traced, see discus- 
sion of binary states in [18]. Another strongly-coupled 
system is produced experimentally, by cold trapped 
atomic gases, by tuning the scattering length to very 
large values via Feshbach resonances. Hydro-type flow 
and small damping of oscillations observed for these sys- 
tems indicate liquid-like properties of matter in this case 
as weU [19]. 

The goal of this paper is to formulate how one can 
tell the difference between the wQGP (with light and 



weakly interacting quarks) from the sQGP picture out- 
lined above, by using a dilepton penetrating probe. More 
specifically, we will discuss the following questions: 

(i) Is it possible to measure the interaction between quark 
quasiparticles directly, using the dilepton probe? 

(ii) Where are the most significant differences between 
sQGP and wQGP scenarios, in respect to the dilepton 
production? 

(iii) How large are they and whether one has a chance to 
observe them, at RHIC and elsewhere? 

We identify two such mass regions: 

(i) M ^ 1.5 — 2 GeV corresponding to T « T(zcro bind- 
ing) and the edge of the quasiparticle continuum, at the 
initial QGP at RHIC. 

(ii) M ~ 0.5 GeV, where we expect to see the contribu- 
tions of the modified vector mesons at T « Tg. 

'\2mt, 
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1 T/T^ 
FIG. 1. Sclicmatic T-dcpendcncc of the masses of qq states. 
A, V, S and PS stand for axial, vector, scalar and pseudoscalar 
states. The dash-dotted line shows a behavior of twice the 
quasiparticle mass. Two black dots indicate places where we 
hope the dilepton signal may be observable. 

The basic idea is very simple: the probability of qq 
process is enhanced by the initial state attractive 
interaction. Attractive interaction obviously correlate q 
and q in space and increases the chances to find qq close to 
each other and annihilate. In such general form, the en- 
hancement persists whether the potential is deep enough 
to make bound levels or not, and whether quark quasi- 
particles have small or large width. 

In the case of pure Coulomb interaction this phe- 
nomenon is well studied. The relevant parameter z = 
ira/v contains the ratio of the coupling strength to ve- 
locity, and the enhancement is in this case given by the 
well known Gamow factor '^'^ 



**We thank V.Zakharov who suggested it to one of us (ES) 
as a possibility, in a discussion after his talk. 



ttThe sign in exponent is the opposite to that on the original 

Gamow factor, for alpha-particle interaction with nuclei or 
that used for HBT correlations, corresponded to a repulsive 
Coulomb force. 
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Note that the result is ^ 1/w at small v (large z) and 1 
at small z. It cancels the velocity in the phase space and 
makes the production rate to jump to a finite value at 
the threshold. 

A good example of the QCD-based color-Coulomb en- 
hancement is the cross section for the production of the 
top quark pair tt, see [20,21]. The enhancement calcu- 
lated in these works is found to be quite significant, in 
spite of the fact that rather large width of the top quarks 
precludes them from forming topponium bound states. 
We will use methods developed in these works in section 
lib. 

In section II we will discuss the dilepton production in 
QCD at T > Tc- But before we get to discuss details, 
let us now point out our main ideas. In Fig.l (which 
is a modified version of a figure present already in [11] 
we show schematic temperature dependence of lowest qq 
states. Dileptons can only come from colorless vector 
pairs (marked by V) . 

New (sQGP) part of this picture is at T > Tc, which 
starts and ends at points marked by small black circles. 
At any given T there are two peaks in the spectral den- 
sity, corresponding to invariant masses of T, L compo- 
nents. On top of that, there is a threshold bump at 
2m* (T), which exists even at T so high that no bound 
states exists. Unfortunately, one can only observe a sig- 
nal integrated over the expanding fireball, or over all tem- 
peratures involved between the initial and freezeout ones. 
The integral tends to wipe out peaks at intermediate lo- 
cations, unless there are special reasons for them to sur- 
vive. Two of the black dots, at the mixed phase T — Tc, 
have a benefit of long time spent there during expan- 
sion. The same is true for the third black point, due 
to near-constancy of the mass of weakly bound states at 
T = 1.5 — 2Tc. The near-threshold bump is at about the 
same location at higher T as well. Thus one may hope 
(and we will show it below) that the structures corre- 
sponding to these endpoints (black circles at Fig.l) may 
be detected. 

Although at finite temperatures and non-zero momen- 
tum relative to heat bath are split into distinct longitudi- 
nal and transverse (L, T) modes, those should coincide at 
zero momentum. Since all masses are large compared to 
relevant T at the time, only pairs with small momentum 
are actually produced. Furthermore, lattice experience 
of similar quasiparticle modes indicate that they follow, 
at least approximately, the vacuum-like dispersion law 
uj^ = -f M^, and so using the invariant mass rather 
than energy would take care of canceling the momentum 
dependence. 

In Fig.l we have shown only states made of u,d 
quarks, ignoring the strange one. Needless to say, 
those exist and can be also produced. The peaks of 
ss ^-like states should be shifted in mass upward by 
2ms ~ 250 — 300 MeV, but their contribution to dilep- 
ton spectra is proportional to square of the electric charge 
q^ = 1/9, which is 3 times smaller than the average of 
q^ = 4/9, = 1/9. Strange states are more promising 



for pseudoscalar/scalar signals to which we turn below. 

II. GENERALITIES 

A. From annihilation rates to the observed dilepton 
spectrum 

Before we proceed to realistic rates, based on lattice- 
based quasiparticle interaction, let tis explain few impor- 
tant points by using a muc;li simpler model. 

As explained above, wc; expect certain structures 
("bumps") to exist in the annihilation rate, related to 
bound states as well as near-threshold enhancement. 
Their exact shape is determined both by the interaction 
between the annihilating quarks (to be studied in detail 
in this work) as well as the interaction with surrounding 
matter. The latter leads to a "width" which we will leave 
to be studied later elsewhere. 

Because of this, it is instructive first to show how sensi- 
tive is the resulting shape of the observable structures to 
such widths. Let us illustrate the point by the following 
example. As a representative of a (near-threshold) bump 
we start with the rate in which the "formfactor" F (2) is 
written as 

F,umr,/2A = Mo ^ ^^^y/,^ + ^ ^^^^j ^ ^ (6) 

In Fig. 2(a) one can see three examples of such bumps, 
with different widths T. 

Although the width is different, the integral over the 
bump is kept the same. Furthermore, this integral is 
normalized to "missing" strength of the spectral density 
due to absent annihilation rate between M = and M = 
Mq (see the second term). We did so because one has 
to be consistent with the so called quark-hadron duality, 
leading to conservation of the total area of the spectral 
density. Provided that one is using a sufficiently simple 
ansatztt, as we do now, the duality restriction provides 
a valuable relation between the bump strength and the 
position of the threshold Mq = 2Mq. 

(Discussion of this and higher order duality relations 
in vacuum correlators can be found in the classic paper 
on the QCD sum rules [22], for more general discussion 
of duality see [23] and for finite T analogs of it as well as 
Weinberg-like sum rules see [24].) 

Here comes the main point we want to make in this 
section. Although three bumps in Fig. 2(a) look very dif- 
ferent, those are for fixed T. In experiment we have to 
integrate the rate over the expanding fireball: let us see 
how these bumps will look after it is done. 



Otherwise just one sum rules does not have a predictive 
power, as experience of QCD sum rules had shown. 
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Although the space-time evolution of the QGP phase 
at RHIC is complicated hydro explosion, here and below 
we will compare with the prediction for RHIC using the 
same parameterization as used by Rapp in [4], concen- 
trating only in the QGP and mixed phase. The physical 
basis for such simplification is that the invariant mass 
of the dilcptons is the same in all frames, and therefore 
motion of the matter can be ignored. (It would not be 
possible if e.g. one wants to calculate the transverse mo- 
mentum of dileptons, or other non- invariant property.) 

We sketch few details here for completeness. The ther- 
mal rate is convoluted with the space time history by 

= / drV{r) / (7) 

The t-dependent volume is modeled as expanding cylin- 
der 

V{t) = (zo + v,t)TT{ro + l^^tY (8) 

where rg — 6.5 fm is the initial transverse overlap, zq — 
0.6 the initial longitudinal lenght, = 1.4c is the relative 
longitudinal expansion velocity of the fireball edges and 
a± = 0.035c^//m is such that at typical freeze-out time 
tf = 15 — 20/m/c. we get a final transverse velocity v± = 
0.6c. The initial temperature of the system is set to T = 
370MeV. The expansion is considered is-entropic with 
the entropy density of the QGP phase given by 

47r2 

SQGP - (16 + 10.5iV/)—T3 (9) 

This fixes the time dependence of the temperature, lead- 
ing to the transition at t ^— Afm/c. The mixed phase 
takes place between Afm/c and 9/m. The fraction of 
QGP in the plasma phase is also calculated from stan- 
dard requirement of constant entropy 

Stot/V{T) = SQGp{T,)fQGP + shg[Tc){1 - fQGp) (10) 

In doing the convolution, one has to define how the 
parameters of the rate, P, Mq depend on T. The results 
shown in Fig. 2(b) correspond to P being independent on 
T while the mass depending linearly on it, Mq = 6T. One 
can see from a comparison of the input to output figures 
that three cases, although still distinct, look much more 
similar. The reason for that is that we kept the integral 
over the spectral density constant, conforming to parton- 
hadron duality. The sensitivity to (unknown) width of 
the bump is reduced because of the averaging. 

The effect of time averaging is maximal in this example 
because we took a simple linear dependence for threshold 
Mq — 6T. In reality, the dependence is nonlinear with 
a minimum (see dash-dotted curve in Fig.l). Although 
the precise shape and position of zero binding point are 
not yet known (better lattice results, please...) one may 
think that real dependence of Mq on T is weaker and the 
effect of averaging be less significant. 




FIG. 2. (a) An example of spectral densities, all of the 
satisfying the duality condition, before time integral, (b) 
This is how the corresponding spectra look like after the 
averaging over the expanding fireball is performed. The lines 
marked by crosses, squares and circles correspond to most 
narrow, the middle and the widest peaks in (a), respectively. 

B. Non-relativistic approach and Green functions 

As mentioned in the Introduction, the main goal of this 
paper is to find out how the strong interaction among 
the initial quarks modifies the dilepton production. Such 
calculations simplify significantly if one can use the non- 
relativistic approach. This is clearly possible for invariant 
masses close to threshold. 

The modification is the same in the inverse process, 
and annihilation of e^e~ pairs into qq has been studied in 
detail in leptonic colliders. The nonrelativistic methods 
we use benefit in particular from discussion of Coulomb 
enhancement of tt pairs near threshold in [20], [21]. 

In [21] the cross section a{e'^e'^tt) is analyzed to 
leading-logarithmic order in QCD in the non relativis- 
tic limit (close to threshold). Their main result, already 
used in [20] , is the modification of the threshold factors 
from the well known leading order expression 
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to 



a = 



3s 



(12) 



Where E is the center of mass energy and Ft is the width 
of the top quark (related to hfetime due to weak decays). 
The relation for dilcpton production rate is obtained by 
obvious substitution of 3G/s into Fq in (2). 

As a further clarification, in the limit F = we can 
rewrite 

9GB+io+ (0, 0) = ^ MO)MO)^S{En - E) (13) 

Note that the sum runs over all states, including scatter- 
ing states, where the sum should be replaced by an inte- 
gral. This expression explains quite clearly the connec- 
tion between the Green's function and the standard non 
relativistic formula for the annihilation of bound states 
(13) gives the flux factor that multiplies the annihilation 
cross section at zero momentum 



a = agg{s = 4mg)|V'n(0)|^ 



(14) 



The non-relativistic Green's function GB-i-irt(?", ?") 
obeys the usual Schrodinger equation: 

[--V^ + V{f} -{E + iT)]GE+ir{r, f) = S^{r- f) (15) 

m 

with inter-particle potential V{r). Analytic expressions 
for Green's function for the pure Coulomb potential are 
well known, see e.g. [20]. For realistic lattice-based poten- 
tial we will use a numerical method following [21], which 
is valid for all potentials less singular than 1/r^ . The ben- 
efit of it is that one avoid summations over levels, and 
bound states are automatically included together with 
scattering ones. As a test, we checked that it reproduces 
well known results for Coulomb potential with accuracy 
at least few per mill. 

The Green's function, thus, contains, at the non rela- 
tivistic level, all the information needed about the states 
of the two particle system. This is important because 
we arc interested not only in bound states but also in 
modification of scattering states, especially close to zero 
energy, to be responsible for the threshold effects we will 
study. The cross section oi q + q ^ + e~ is calculated 
by simply inverting the previous one from phase space 
considerations (neglecting the lepton masses) 



cr, 



.(s) 



1 Cri++/ ~q^q{s) 



(1 



Ami 



) 



(16) 



to solve the Schrodinger (or Klein Gordon) equation for 
different temperatures. 

Another ingredient taken from lattice calculation is the 
quasiparicles masses. In [16] Petreczky et al. have shown 
that the masses of the quarks in the region of 1.5 — 3Tc 
stay roughly constant at a value of Mq « IGeV As in 
[12] we fix this value all the way up to Tc, being aware 
that the quasiparticles should become heavier close to 
the transition temperature. 

It was also pointed out in [12] that in order to trig- 
ger the phase transition at T = T^. the lowest qq bound 
state is required to become massless. The previous static 
potential is not enough to achieve this, and a quasi- 
local interaction induced y the instanton-antiinstanton 
molecules [12] [13] is needed. We will model this interac- 
tion by the following potential: 



^inst — 



(r2 + p2)3 



(17) 



where p = 1/3/to is the typical instanton size. The value 
Uq is fixed at by the requirement of the lowest level to 
be exactly massless (in the chiral limit). At higher tem- 
perature, the instantons are suppressed, which we model 
with the following damping factor 



n(p,T) =n(p,Te)exp(-i(27Ve 



+ Nf){7rpnT^-T\)) (18) 



Finally, the polarization of the instanton molecules along 
the (Euclidean) time direction leads to a distinction of 
the TT and p channels, see discussion in [12]. To take 
this into account, the authors of [12] have multiplied the 
strength of the non-local interaction by a correction fac- 
tor Uq — > UoF where for vectors F = 0.7. 



III. RESULTS 

First, let us explain our units. Since the effective mass 

of quark quasiparticle in the temperature interval con- 
sidered is not known accurately, we use it as our basic 
energy unit. In plots twice this value appear as a thresh- 
old, to which we tentatively ascribe to it a particularly 
simple value 2Mq — 2 GeV: the reader should however 
be warned that this is a "GeV" in quotation marks, to 
be rescaled appropriately later when the value of quark 
effective masses in QGP will be better known. 



A. A near-threshold region 



C. Lattice- based potentials for static quarks 

As in [13] we use a potential extracted by parameter- 
izing the Bielefeld data [15] on the effective static po- 
tentials (for details see [13]). These potentials are used 



Although the underlying potentials describing interac- 
tion between quasiparticles in the channels we consider 
arc always attractive, due to presence of bound/virtual 
states the effective scattering amplitude is complex and 
its real part may even change sign. In particular, it is well 
known from quantum mechanics that since a scattering 
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amplitude changes sign at energy equal to the position 
of the level, an overall attractive potential can cause ef- 
fective interaction of positive energy particles to be re- 
pulsive, if the level is close to zero, However this effect 
holds only for scattering states, or particles moved far 
away from the origin (where the potential is negligible). 
We however are interested in a shape of the annihilation 
signal: and thus we are only concerned about the wave 
function at the origin. Even in a case when a level ap- 
proach zero, the effect of the attractive potential is only 
to increase the annihilation rate: there cannot be any 
opposite effect. 

All of it can be studied in the simplest problem possi- 
ble, that of a spherical potential well: 



V = 



0, (r > a) 
-Vq, (r < a) 



(19) 



we find levels at thresh- 



Parametrizing Vq - 4^2^ 
old {E = 0) for n integer and n is the number of levels 
below threshold. It is straightforward now to solve the 
Green's function equation for this potential (15). Again, 
as we are only interested in the Green's function at the 
origin, we only need to calculate the s wave equation 
with the standard boundary conditions. The result of 
this simple calculation (in the limit F — > 0+ is 



3G 



E+M+ 



(0,0) 



1 



4 7r((a;2 - fc2) cos'^{uja) + k^) 



(20) 



And = M{E + VJ,), = ME. When n is not 
an integer, the Green's function close to zero behaves as 
— 4 — , that is, the value of the free Green's function. 
When n is integer, we have a level at i? = what leads 
to a divergence. We see then that the modifications at 
threshold come from the presence of bound states and 
virtual levels close to zero binding. In the case of coulomb 
potential, the finite value of this object at threshold is a 
consequence of the infinite number of states close to that 
energy. 

We also observe that, according to (12), the cross sec- 
tion at any energy is bigger (or equal) than the free one. 
as correspond to an attractive potential. It is interesting 
to note that the oscillations at positive energy correspond 
to virtual levels^ as the continuation of the denominator 
of 12 to imaginary values of k is the condition for bound 
states in the potential. The position of the maxima, then, 
reflects the position of these virtual levels. 

Repeating the calculation that lead to (4) (see, for ex- 
ample, [5]) with the cross sections expressed in terms of 
the Green's function (12) we obtain the following expres- 
sion 



9ne™(M) = -43GM-2™,+zr,(0,0) 



(21) 



The calculation of the modification of the spectral den- 
sity is now straight forward. The parametrization of lat- 
tice results plus the non local interaction described in the 



previous section gives us a close expression for the inter- 
particle potentials at different temperatures that we can 
use to solve (15). We do this numerically by following the 
numerical method of [21] based on finding two indepen- 
dent solutions of Schreodinger equations. We evaluate 
SGM-2m,-i-ir,(0,0) directly for "realistic" lattice-based 
potentials. 




FIG. 3. S>G£;+io+ (0> 0) versus energy/mass for a=l also in 
inverse mass units: Free case (solid), n=0.01 (cross), n=0.5 
(box), n=0.9 (circle) 

We still have to address the issue of the width. It is 
clear that the quasiparticles will have a thermal width in 
the equilibrated medium. Perturbative calculations for 
this quantity exist in the literature and to lowest order 
they are (at rest) 



T = a{N,Nf) 



4-77 



(22) 



With a{N,Nf) « 1.40 for iV = 3 and Nf = 2. For 
T « (l-3)Tc and a,, « 0.5 gives a width T « OA-lGeV. 
These values of the width, comparable to the quasiparti- 
cle mass, are too big and would make very complicated 
to be able to extract a value for the masses from lat- 
tice data. Of course, there is no reason to believe than 
the perturbative expression makes sense for such a big 
value of ag. That is why we follow a phenomenologi- 
cal approach, leaving the width as a free parameter and 
presenting our calculations for different values. 
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FIG. 4. Modification of tlie spectral density at Tc witfi re- 
spect to invariant mass in units of Mq. The peak at 1600 
MeV corresponds to the p (non relativistic approximation). 
The peak at 2000 MeV reflects the threshold enhancement 

We now move to present the results, starting with the 
worning. In figure 4 we show the modification of the 
spectral function, (the ratio (21) over (3)) at T = Tc 
(r — IQMeV) . One can see the threshold enhancement 
at 2 GeV and the peak around 1600 MeV corresponding 
to the bound state that should represent p aXT = Tc- Al- 
though the potential we used has the same suppression 
of instanton molecules as in [12], instead of relativistic 
Klein Gordon equation used in that work we now use 
the nonrelativistic Schrodinger equation and thus found 
much less deeply bound state, with only about 400 MeV 
binding. (We checked that Klein Gordon equation with 
the same potential modified by for Schrodinger and by 
simulating the velocity-velocity term by effectively dou- 
bling the couphng [12] [13] indeed shifts p much lower, to 
about 600Afey, as found in that paper.) 

Thus at r « Tc the relativistic effects are very impor- 
tant. It means that one should either (i) develop fully 
relativistic theory of Green functions and dilepton rates, 
or (ii) restrict the discussion to the nonrelativistic do- 
main at somewhat higher T not too close to Tc- In this 
work we only follow the latter option: the reader thus 
should be aware that our results should not be trusted 
close to Tc- 

What we learned from these calculations is how the 
spectral density changes as the bound states become less 
and less bound until the system arrives to zero binding 
(zero binding point). In figure 5 we show the shape of 
the dilepton spectral density for different temperature 
(1.2 Tc, 1.4 Tc, 1.7 Tc, 3 T^). One can observe how 
exactly enhancement in the bound state and threshold 
region changes. The hight of this enhancement depends 
on the proximity of the level; at 1.7 Tc, where the bound 
level is close to threshold, we observe a big modification 
of the spectral function at 2GeV. Note that proximity of 
the level to threshold happens in rather wide T interval, 
roughly between 1.5 Tc to 2 T^ (the zero binding point 
[11]. Note also that the enhancement is still seen at tem- 



peratures as high as 3 Tc, when all bound states have 
already been melted. 

We show in figure (6) the dependence of the modifica- 
tion factor with the width of the quasiparticles at a fixed 
temperature of 2 Tc. The general features of the modi- 
fication are the same. However, the extra enhancement 
due to the proximity to threshold is very sensitive to the 
width of the quasiparticles. But, we expect at least a 
factor 2 from the final calculation of the rates with the 
standard expression (3) for any reasonable value of the 
width. 



B. Dilepton signal at RHIC 

In this section we discuss the time-integrated dilepton 
production rates. We start with the QGP phase. We can 
observe an enhancement in the region close to 2 GeV. 
This is because in the whole region of temperatures ac- 
cessible at RHIC we observe the bump in the cross section 
due to threshold effect. The contribution below thresh- 
old is strongly suppressed due to the mass effect. The 
value observed comes from the bound state moving from 
threshold till the value at Tc. 

When we look at the contribution from the mixed 
phase, we can observe a significant threshold contribu- 
tion. However the main contribution is given by the 
bound state at Tc. We remark again that this contri- 
bution should not be considered too seriously due to the 
relativistic character of the bound states. 
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FIG. 5. Modification of the spectral density versus the in- 
variant mass in Mq units for different temperatures: (a) 1.2 
Tc (cross) and 1.4 Tc (circle), (b) 1.7 Tc (cross) and 3 Tc 
(circle) and the correction due to quark mass (line) . 



t 




FIG. 6. Modification of the spectral density versus the 
invariant mass in Mq units . Dependence on the width 
(at T = 2Tc)\ r = lOMeV (cross), F = 50MeV (box), 
F = IQQMeV (circle) 

Finally we show the total modification of the produc- 
tion due to the QGP. We show the effect of the widths. 



As the contribution of the mixed phase is much smaller 
than the QGP one, the effect of the bound states is not 
as prominent. At the same time the contribution from 
the threshold produces an enhancement that survives 
the time integration. Although both effects depend on 
the assume widths, they are present for values as big as 
F — IQQMeV . Thus, we do expect an observable modi- 
fication of the spectrum in sQGP, relative to the wQGP 
calculation by Rapp [4]. 




FIG. 7. The contribution of (a) the QGP phase and (b) 
the mixed phase at RHIC to total dilepton yield, normalized 
to the "standard rates" (note that both approach 1 at high 
masses). The calculation is done for F = IQMeV), the invari- 
ant mass is given in units of Mq. 

IV. TWO-PHOTON SIGNALS 

At first glance, 77 spectroscopy looks even more 
promising than dileptons, since such decays can hap- 
pen for wider range of channels - such as pseudoscalars, 
scalars and tensors. 
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1.6 1.8 2 2.2 2.4 2.6 2.8 

FIG. 8. The ratio of dilepton production from the 
QGP+mixed phase to "standarad rate" , vs the invariant mass 
(in Mq units). The effect is shown for two widths (per quark), 
(a) r = lOMeV, (b)r = lOOMeV 

Pseudoscalars: In vacuum, as well as the hadronic 
phase T < Tc the pseudoscalars tt, are Goldstone 
bosons, and their masses are protected by chiral sym- 
metry, and thus if all quarks be massless they should 
remain zero. As it is wel known, the near-SU(3)-singlet 
•q' is not a Goldstone boson due to chiral anomaly and 
instantons, which explicitly violates the chiral U{\) sym- 
metry. A sketch of levels motion with T is shown in Fig. 9. 
As T grows, the mixing between 77' and t] is expected to 
change as well, they start at T = close to pure SU(3) 
singlet and octet, but above Tc they eventually rotate 
into purely light wu, dd and strange ss states. This is 
related with the fact that the role of instanton-induced 
interaction deminishes with growing T, both because of 
reduced wave function at the origin of the bound states 
and because of instanton suppression by Debye screening. 

For estimates of the rate we will need the correspond- 
ing electromagnetic widths, given by their total widths 
and diphoton branching ratios from the Particle Data 
Tables: 

T^^^^ = l.l?>keV * .m= AQkeV (23) 
r,,-_^^ = .2 MeV * 0.02 = 4 keV (24) 



Scalars: As it is well known, the situation with scalar 
mesons is rather complex. The debates of whether the 
lightest scalars such as /o(600) (known formely as a 
of Gel-Mann-Levy sigma model) is a qq state go back 
decades, and some heavier scalars are now believed to 
be partly diquark-antidiquark states, and also mix with 
the scalar glueball. Not going into this discussion we only 
comment that above Tc (in the chiral limit) there must be 
degeneracy between pseudoscalar and scalar chiral part- 
ners of the pion and eta mesons. In other words, if quark 
chirality is conserved, there are identical states made of 
qLqB. and qjiqL quarks. 

Tensors: /2(1270) is a famous example of a state seen 
in 77 collisions, its partial width is comparable to that 
of scalars and pseudoscalars, namely 

r/_^^ = 185 MeV * 1.410"^ = 2.6 keV (25) 

However aiT > Tc the situation is radically different: the 
calculations made in [13] show that all^^ P-wave states 
seem to be dissolving very close to Tc, and in fact the 
lattice-based effective potentials lead to a single s-wave 
bound states. 

An estimate of the production rate may be done as 
follows 

A^^^ = j T^^ndtdV (26) 

where the integral is done over space-time occupied by 
sQGP. Using for estimate tqcp ^ h fm ^ 1/(40 MeT^) 
times gives small parameter TIqcp ~ (4 keV/ 4QMeV) ~ 
lO"**, originating ultimately from the electromahnetic 
coupling ttg^. The integral of density over volume ndV 
is dominated by later stages and thus is approximately 
equal to the particle number produced at chemical freeze- 
out. Thus we conclude that for 77,77' 2-photon signal the 
in-matter (modified) signals are about 4 orders of mag- 
nitude below those coming after freezeout. 

If instead of T w Tc we think now about the contribu- 
tion to the time integral of the early times, at the vicinity 
of the zero binding point, the temperature is about twice 
higher but the mass of the state we would like to produce, 
with the energy w 2Mq k, 2 GeV is also twice higher than 
the mass of the 77', resulting in about the same Boltzmann 
factor. 

Good news about these rates is that there should be 
more states similar to tt, 77: adding scalars and states 
containing "plasminos" instead of quark quasiparticles 
leads to additional factor 2*4=8 in the rate. 

Unfortunately there are also small parameters in- 
volved, leading to further troubles for the 77 signal. The 



^^The only exception is the scalar glueball channel in which 
color Casimir is factor 9/4 stronger and P-wave state sur- 
vives for a while. However glueball states are not relevant for 
electromagnetic probes anyway. 
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widths of the pseudoscalar mesons in vacuum are 
proportional to etc which (as one of the order pa- 

rameters of chiral symmetry) must vanish at T > Tc (in 
the chiral limit). What it means is simply the statement 
that those particles are bound state of LR or RL (where 
we mean left and right handedness of quarks) and in the 
chirally symmetric QGP phase L and R are like differ- 
ent quark flavors, which simply cannot annihilate each 
other as they are not antiparticle of each other. In order 
that to become possible, one should admit chirality flip 
in the process, or the first power of a current quark in 
the amplitude, resulting in additional suppression factor 

r^^(T > Tc) . rris 2 leyj\ 

which reduces the effect by about an order of magnitude 
for strange quarks (e.g. in rj.rj') and much more so for 
the pion. 

Although the authors are not qualified to tell what 
exactly the limits of experimental detection of this effect 
is, it looks to be rather difficult to measure these effects. 
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FIG. 9. Schematic T-dependence of the masses of qq pseu- 
doscalar states. Two dash-dotted lines shows a behavior 
of twice the quark quasiparticle mass, for light and strange 
quarks, respectively. Two black dots indicate places where 77 
signal may be observable. 

V. SUMMARY AND DISCUSSION 

We have proposed in this work new direction of dilep- 
ton experiments, namely looking for the bound vector qq 
states and threshold enhancement in the QGP domain. 
We in particularly identified 2 regions of the dilepton in- 
variant mass as most promising, corresponding to two 
endspoins of the sQGP domain. 

In a nonrelativistic approach, we worked out a numer- 
ical method for the determination of the Green function 
and annihilation rate, using all states, bound and un- 
bound. We used realistic interaction and have studied 
in detail the shape of spectral density as a bound level 
approach a threshold. 

We have also performed realistic averaging over an ex- 
pansing fireball, using parameters for RHIC, and demon- 



strate that this evraging does not eliminate the contribu- 
tion of the peaks near two ends of the sQGP region. The 
actual observability of states still depend on the issue of 
the width of these states. We have decided at this point 
to leave this issue open for further studies. The reason 
for that is the rather uncertain theory of sQGP and its 
near-perfect liquid properties. For the first investigations 
of the "photo-effect" -like reactions with splitting of the 
binary bound states see [25]. Lattice studies do see these 
peaks, which is encouraging, but the maximal entropy 
method they use in [14] was not studied enough to tell 
how accurately the widths can really be estimated by it. 

We pointed out that in principle a dilepton study can 
be complemented by a 77 one, for scalars and pseu- 
doscalars from sQGP. Interesting signals can be expected 
around Tc and the same threshold phenomena near the 
endpoint as for dileptons. However, much higher back- 
ground makes it more difficult to do. 
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